Vortex flow for a holographic superconductor 
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We investigate energy dissipation associated with the motion of the scalar condensate in a holo- 
graphic superconductor model constructed from the charged scalar field coupled to the Maxwell 
field. Upon application of constant magnetic and electric fields, we analytically construct the vortex 
flow solution, and find the vortex flow resistance near the second-order phase transition where the 
scalar condensate begins. The characteristic feature of the non-equilibrium state agrees with the 
one predicted by the time-dependent Ginzburg-Landau (TDGL) theory. We evaluate the kinetic 
coefRcient in the TDGL equation along the line of the second-order phase transition. At zero mag- 
netic field, the other coefficients in the TDGL equation are also evaluated just below the critical 
temperature. 

PACS numbers: 11.25.Tq, 74.20.-z, 74.25.Qt 



I. INTRODUCTION 

Much attention has been given to the application of 
the AdS/CFT (anti-de Sitter/conformal field theory) du- 
ality [H to condensed matter physics after discovery of 
holographic superconductor models 0, Since the 

AdS/CFT duality is a valuable tool for investigating 
strongly coupled gauge theories, the application might 
offer new insight into the investigation of strongly in- 
teracting condensed matter systems where perturbative 
methods are no longer available. 

The holographic superconductor model constructed by 
charged scalar condensate 3 is classified into type II su- 
perconductors, as it possesses vortex solutions 0-01 in 
a background magnetic field. Furthermore, it has been 
shown that a triangular vortex lattice solution is the 
most favorable solution thermodynamically just below 
the second order phase transition at long wavelengths 
0. As already seen in Refs. [sl-fioj. these equilibrium 
states are described by the Ginzburg-Landau (GL) the- 
ory. This suggests that non-equilibrium states of the 
holographic superconductor in the background magnetic 
field are also described by the time dependent Ginzburg- 
Landau (TDGL) theory. Indeed, it has been observed 
that the dynamics in the absence of magnetic field is de- 
scribed by the TDGL theory [iMl- 

Motivated by this, we investigate the non-equilibrium 
steady state of the vortex lattice solution 0] in the pres- 
ence of a small constant electric field E. According to 
the TDGL theory, the vortex flows at a constant velocity 
in a direction perpendicular to both the magnetic and 
electric fields so that the Lorentz force on the vortex is 
balanced by the background electric force. The energy 



'Electronic address: 
t Electronic address: 



maeda302Qsic.shibaura-it.ac.jp] 



tokamura® kwansei . ac . j p l " 



dissipation associated with the vortex motion occurs in 
the core of the vortex (vortex flow resistance), as the 
superconducting state disappears there. 

In the TDGL equation, the evaluation of the kinetic 
coefficient F (for example, see Ea. (|B3[) ') is important 
in observing the dissipation process or the spectrum of 
quasi-particlcs around the core of the vortex. While it 
is generically difficult to derive the coefficient from the 
microscopic point of view in strongly interacting con- 
densed matter systems, it can be evaluated in the holo- 
graphic model. So, it is interesting to explore the dissi- 
pation mechanism associated with the vortex motion in 
the framework of the AdS/CFT duality. 

In this article, we perturbatively construct the vortex 
flow solution as a series expansion of the small electric 
field E and derive the R-current just below the criti- 
cal temperature where the scalar condensate begins. We 
find that vortices flow at a constant velocity and that 
the Ohmic dissipation occurs associated with the vortex 
motion, as predicted by the TDGL theory. The kinetic 
coefficient F is evaluated along the line of the second or- 
der phase transition in comparison with the TDGL the- 
ory. In the absence of magnetic field, we also derive the 
other coefficients in the TDGL equation from the value 
of the the scalar condensate, the correlation length, and 
the London equation Q. 

The plan of our paper is as follows: In Sec. II, we ex- 
pand equation of motion for the scalar field as a series 
in E. In Sec. Ill, we perturbatively construct the vor- 
tex flow solution by Green function method. In Sec. IV, 
we derive the net R-current by solving Maxwell equation 
and evaluate the kinetic coefficient. We briefly review 
the TDGL theory in Appendix B, and the other coeffi- 
cients in the TDGL equation are derived in Appendix C. 
Conclusions and discussion are devoted to Sec. V. 
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II. BASIC EQUATIONS IN 
EDDINGTON-FINKELSTEIN FORM 

We consider the (2+l)-dimensional holographic super- 
conductor model described by a dual gravitational theory 
in four dimensions (AdS4) coupled to a charged complex 
scalar field ^ and a Maxwell field Afj, [3, 0] ■ For simplic- 
ity, we take a probe limit where the backreaction of the 
matter field onto the geometry can be ignored 0] . 

The background metric is given by AdS4-Schwarzschild 
black hole with metric 



{~h{u)dt' + dx^ + dy") . 

u a[u) 



h{u) = l-u^ 



a[T) = 47rr/3 



(2. fa) 
(2.fb) 



where L and T are the AdS radius and the Hawking tem- 
perature, respectively. We take the coordinate u such 
that the AdS boundary is located at u = and the hori- 
zon is set to be M = f . 

Under the probe limit, the action of the matter system 
S = {L^/2Kle'^)S is written by 



p-2 



ml* I 



(2.2) 



where m and e are the mass and charge of the scalar field 
respectively, and 



(2.3) 



Hereafter, we consider the action (|2.2p in the simple case 
rn?L^ ~ —2. The equations of motion are given by 



L2 



f = 



(2.4a) 
(2.4b) 



For a gauge choice, we choose a gauge Au = in the 
metric (j2.f p . The asymptotic behavior of 4* and (/i ~ 
t, X, y) near the AdS boundary are 



^' ~ ci{t, x,y)u + C2{t, X, y) 



A„ 



Af^{t,x,y) + J^{t,x,y)u . 



(2.5a) 
(2.5b) 



According to the AdS/CFT dictionary, the expectation 
values of the dual scalar operator with conformal 
dimension two and the R-current {fi ~ t, x, y) are 
represented by the coefficient C2 and JT^ as 



O2 



J" 



^/2e SS 
La Scl 
SS 



V2L^ «2 



2K|e 



■ C2 



SAu 



u=0 



2K|e2 



(2.6a) 
(2.6b) 



^ Our definition of ( O2 ) is the same as Ref. [Tsj in units where 
L2/2«;2 = 1. 



respectively. We consider the condensation of the scalar 
operator O2, and impose an asymptotic boundary con- 
dition ci = to eliminate the source term in the dual 
theory. 

Since we are interested in the superconducting region 
just below the second order phase transition, the ampli- 
tude of the scalar field j^*! is very small. So, one can 
expand ^' and A^ in powers of a small parameter e as 

* = ei/2 ^, + o(g3/2) ^ A^ = A^ + ea^ + 0{e^) . 
Then, Eq. (j2.4b|) at 0(e°) is reduced to 



y.F^" = 0, 

where F^^ = 9^^^ - By Eqs. 

of motion for and become 



(2.7) 

, the equations 



(2.8a) 
(2.8b) 



where and /^i, are defined by := ~ i and 
f^j,y := 2d[^a„], respectively. 

We consider a zeroth order solution of Eq. p.7p gener- 
ating a constant electric field E and the (upper) critical 
magnetic field Bc2 at the second order phase transition 
. This is given by the following form: 

At = nil -u), A,^-E{t-u,)-B,2y , (2.9) 

where := J" du/ah{u) and /i is the chemical poten- 
tial. The boundary conditions at the horizon u = 1 arc 
determined by the regularity condition, i.e., At{x, y,u ~ 
1) = and \F'^(x,y,u ~ 1)\ < 00. Our strategy is to 
solve Eqs. p.8p perturbatively for small E under the ex- 
ternal gauge field (|2.9p . In the E = case, Eq. p.8a|) is 
solved as the Landau problem [13) , and Eq. (j2.8bp is also 
formally solved 0. 

Let us expand Eqs. ()2.8p in powers of E as 

(2.10a) 
(2.10b) 



tp = tpo + E t/ji + ■ ■ ■ 
a, = a^)+Ea!^^ + 



As shown later, it is convenient to adopt an advanced 
null coordinate v := t — and a coordinate y defined by 



E 



(2.11) 



Under the coordinate transformation it, u, x, y) 1— >■ 
{v,u,x,y), the metric p.ip and the gauge field ()2.9p are 
transformed as 



ds^ 



h(u) dv dv du + dx 

a 



2E 

Bc2 

Ay — /z(l — m), A. 



dydvj +0{E'^) , 
4 



a hiu) ■ 



(2.12a) 

-Bc2y- 
(2.12b) 
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Thus, the external electric field E appears only via the 
metric in the Eddington-Finkelstein form (|2.12ap . 

If the holographic superconductor obeys conventional 
type II superconductors, the vortex flows at a constant 
velocity —E/Bc2 along the y-direction by the Lorentz 
force. This implies that there is a "static" vortex solution 
in the metric (|2.12ap . So, we take an ansatz for the scalar 
field V^: 



ip{x, y, u) = 



dp 



^iy,u;p) = ^o{y,u;p) + E^i{y,u;p) -] . (2.13a) 

Defining the differential operators £p and C as 



1 



9^ 

dy^ 



{p + B,2y? 



C 



d h[u) d 



du du h{u) 



(2.14a) 
4 ' (2-14b) 



we obtain the equations of motion for ^nd from 
Eq. (|2?8a| as 



-Cp Co — , 
2 d 



£p6 = j , 
d 1 i/i(l — u) 



Co 



aBci dy \ du u ah{u) 
The boundary conditions of ^ are represented by 

C2{x,y)u^ (w ^ 0) 



(2.15a) 
(2.15b) 



^{x,y,u) = 



regular 



(2.16) 



III. THE CONSTRUCTION OF THE VORTEX 
FLOW SOLUTION 

In this section, we construct the solutions a-nd d of 
Eqs. (|2.15ap . Following the ansatz in Ref. 0, we separate 
the variable as = Pn{u)Dn{Yp), where Yp is defined 

by 



Yp 2Bc2 



B. 



c2 



(3.1) 



The equations for pn{u) and Dn{Yp) arc derived from 
Eq. (|2.15ap as 



l-'Pn\u) = — pn[u) 



(3.2a) 



^ In Ref. Q , the expansion along the x-direction is given by Fourier 
series rather than the Fourier transform. As shown in the Ap- 
pendix A, however, we can formally write the vortex lattice so- 
lution in the framework of the Fourier transform. 



where A„ is a separation constant. The solution of 
the equation (j3.2b[) satisfying the boundary condition 
lim|y|^oo \Dn{Y)\ < oo is given by 



A„ = 2?i + 1 



(71 = 0,1, 2,...) , 



1 



1/2 



Y 



/27r2"n!y " \ J2 . 
where Hn is the Hermite function defined by 



(3.3a) 
(3.3b) 

(3.4) 



The function Dn{Yp) in Eq. (j3.3bp is the rt-th en- 
ergy eigenfunction of a harmonic oscillator centered at 
y = —p/Bc2 and it exponentially decays for large \Yp\. 
As discussed in Ref. Q, the upper critical value Bc2 
is determined by n = and the solution po satisfy- 
ing the two boundary conditions (|2.16p was numeri- 
cally obtained. Therefore, we shall adopt ?i = solu- 
tion, Co = Po{u) Df){Yp) as the leading order solution of 
Eq. (pfT5a|) . 

We derive the next order solution Ci of Eq. (|2.15ap 
by constructing Green function. In general, Ci includes a 
component proportional to Co- Hereafter, we shall remove 
this component from Ci because it can be absorbed into 
the leading order solution Co- 

Introducing the inner product for £)„ as 



(C,77) := / dYe{Y)il[Y) , 



{Dn} forms a complete orthonormal set 



J2Dn{Yp)Di{Y;) = 6{Yp-Y;) 



n=0 



S{y-y') 
V2B. 



c2 



(3.5) 

(3.6a) 
(3.6b) 



As well known, {Dn} satisfies the relation 



VtT+T m ~ n + 1 
{Dn, YDm) = { i/n m = n - 1 

\m -n\^l 



(3.7) 



In terms of the complete orthonormal set {Dn}, let us 
construct the Green function Gp(u, y\u',y') of the oper- 



ator £,p in the form 



Gp{u,y\u',y') 



V2B, 



c2 



Y,9n{u,u')Dn{Yp)Di{Y;) . (3., 



n=l 



Suppose that the two point function g„ (n = 1, 2, - - - ) is 
the solution of the equation 



^£ gn{u,u' ) ^ S{u ~ u') (3. 



9) 
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satisfying the boundary condition (|2.16p . Then, using 
the completeness of £'„(F) in Eq. (j3.6b[) . we find that 
Gp satisfies 

£pGp{u,y\u',y')^S{y-y')S{u-u') 

- y^2S{u - u' ) Do{Yp) dI{Y;) . (3.10) 

This indicates that Gp is the Green function in the solu- 
tion space orthogonal to the component on Dq. 

In terms of the Green function Gp, the formal solution 
of ^1 is represented as 

^i{y,u]p)= I du'dy' Gp{u,y\u',y')j{y',u';p) . (3.11) 



Substituting Eq. ()3.8p into the solution p.ll|) and using 
the orthogonality condition p.6a|) . one obtains 



f du' gi{u,u') (u'-^ 



■■Di{Yp) 

I d 1 *Ai(l — u') 



Di{Yp) 



ah{u') 

{pR{u)+ipj{u)} , (3.12) 



where and pi are the real and imaginary parts of the 
u'-integral, respectively. 



IV. ENERGY DISSIPATION ASSOCIATED 
WITH THE VORTEX FLOW 

In this section, we investigate energy dissipation 
caused by the vortex flow solution constructed in the 
previous section. As shown below, the R-current asso- 
ciated with the vortex flow agrees with the one predicted 
by the TDGL theory (sec. Appendix B). So, we can eval- 
uate the kinetic coefficient F in the TDGL equation (|B3[) 
from the R-current. The other coefficients are also evalu- 
ated from our earlier results in Appendix C. We begin 
by calculating the R-current induced by the vortex flow. 



R-current 



( ) in Eq. (j2.6b[) can be expanded as a series in e 
near the second order phase transition as 



2K|e2 



(5( 



(4.1) 



where 5 { ) — 0{e) = 0{"i'^). The flrst term is caused 
by the normal fluid, which is independent of the vortex 
motion. Hereafter, we will calculate the subleading term 
(5 ( ) , as it is induced by the vortex motion. 

For simplicity, we shall focus attention on calculating 
the net (total) current 5 { ). We define the net value 
A of a quantity A in the original coordinate (x, y) as 



A := dxdy A. 



(4.2) 



By Eqs. (j2.5bp and ()2.6bp . 5 { J"" ) becomes 



L^ae 
2^ 



(4.3) 



We first derive Qx by solving Eq. (j2.8b| . As in 
Eq. p.lOp . the net bulk current jx can be expanded as 



(4.4) 



The leading term ji^"^ is clearly zero because the current 
circulates (see, Ref. Q) as 



.•(0) 

Jx 



-<^dy\lp0\ 



(4.5) 



(1) 



This implies that the leading term ax in Eq. (|2.10bp is 
zero, and that is written hy ax — E ax 

By Eqs. (|XBal) . (|X7|) . and (|XT^ . the subleading term 
ji^'' is calculated as 



Jx 



(1) 



-2e5 [{Dxi^o)^i + {Dx^pi)Ho] 



2eC 



C :-- 



dp |C(p)p , 



Po{u)pRiu) , 



(4.6a) 
(4.6b) 



which is independent of the coordinate v. In the spirit 
of finding the solution of Eq. (j2.8bp . we assume that aj^^' 
is also u- independent. Then, averaging the x-component 
of Eq. (|2.8b|) and using Eq. we find 



d_ / da^^\u) 
du \ du 



ay/2Bc2u'- 



■Po{u)pr{u). (4.7) 



Here, we used the fact that the spatially averaged value 
in Eq. ()4.2p is zero for the derivative terms of with 
respect to the spatial coordinates, x, y, i. e. , dxa^ = 

dyCL^ = 0. The regularity of ai^^ at the horizon deter- 
mines the solution as 



_a(i)(u)^- 
du ^ ay/2Bc2h{u) 7m 



Po{u')pr{u') 



(4.8) 



Substituting Eq. (|4.8p into Eq. (|4.3p . the net R-current 
becomes 



6( 



eEL^C 



{Po, Pb), 



(4.9) 



where (po, Pr) is the inner product defined by 
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Next, we express the net R-current in Eq. (j4.9|) in terms 
of the expectation value of the scalar operator ( ©2 )• 
Under the boundary condition (j2.16|) . the operator £ is 
clearly Hermitian for the inner product (|4.10[) . So, using 
Eqs. (|3^ . ((3J)) . and ((3?T2|) . we obtain the following 
equality: 



(Po, C-Pr) = (po, pr) + 2{ po, u 

= {LpQ, PR) = — {po, pb) . 



A. El 

du u 



(4.11) 



The boundary condition (|2.16p simplifies the equality as 



(po, PR) = ^2B~^°^^ " ^ 



c2 



(4.12) 



Substituting Eq. (|4J2l) into Eq. (gS]), S { ) is ex- 
pressed by the expectation value of the dual scalar oper- 
ator ( ©2 ) : 



eL^a^ E Cpg(l) _ iW 



2K|e2 B. 



E 

c2 



(4.13) 



Here, 4*0 e^/^ t/jq and the coefficient /3 is defined by 











|C2 


2 



lim„^o |po(w)/it^P 



(4.14) 



Eq. (|4.13|) shows that a finite DC-current is induced by 
the motion of the scalar field in a parallel direction with 
the applied electric field E. Thus, the vortex flow re- 
sistance appears by the vortex motion in the holographic 
superconductor model. In the bulk side, the energy dissi- 
pation (Ohmic dissipation) associated with the resistance 
is represented by the energy absorption of the scalar field 
by the black hole. As shown in Eq. (|D5[) . the energy flows 
into the bulk from the boundary via the external electric 
fleld. It is transformed into the energy of the scalar fleld 
in the bulk. Since the scalar fleld falls into the black 
hole horizon, as it moves in the y-direction, the energy is 
absorbed into the black hole. 




FIG. 1: T-dependence of Z := kIV/{L^ii^) with 11 fixed. 
Here Tc is the critical temperature for 3^2 ~ 0. 



gauge Ay^ = 0, there is still a residual gauge transforma- 
tion 

A^{t, X, y, u) A^{t, X, y, u) + d^A{t, x, y) , (4.16a) 
^{t,x,y,u) e'^'-^'^''^^ ^{t,x,y,u). (4.16b) 



Then, Eq. ()4.16p acts on the source Af^ of the R-currcnt 
and on the condensate ( O2 ) dual to as 

A^,{t, X, y) Af,{t, x, y) + df,A{t, x, y) , (4.17a) 
( 02{t, x,y))^ e'A(*'-'«' ( 02(i, X, y) ) . (4.17b) 

This is a "background local U(l)" transformation of the 
dual field theory, indicating that the gauge coupling con- 
stant e« of ( 02 ) is unity, i.e., e* = 1. 

By solving numerically Eq. p.2ap with n = 0, wc ob- 
tain Z := (a//i/3)^ = K|r/(L^^^), which is a function of 
a/p, (X T/Tc only. In Fig. 1, we present r(T) as a func- 
tion of T. It increases as T decreases from the critical 
temperature Tc at zero magnetic field, Bc2 =0. In the 
Bc2 — ^ limit [T ^ Tc), both a/ p and Z approach criti- 
cal values ad p ~ 0.25 and Z ~ 0.54, respectively, which 
are independent of the critical temperature Tc. Thus, we 
finally obtain F as 



B. The kinetic coefficient F 



The form of the expectation value (|4.13p agrees with 
the averaged value of the current in TDGL theory (jB17p . 
Then, the kinetic coefficient F is given by 



F(T) L^p^ [a/p) 



/32 



2, ,2 



L'p 



(4.15) 



T,B=Ba2iT) 



We can easily show that e, = 1 in the following argu- 
ment. As seen in Eq. (|2.3p . the gauge coupling between 
Af^ and is given in the form, (9^ — z A^)*?. Under the 



in the T ^ Tr limit. 



303.2 X t2 



(4.18) 



V. CONCLUSIONS AND DISCUSSION 

We have investigated the vortex motion of a holo- 
graphic superconductor constructed by a gravitational 
model of complex scalar field coupled to the C/(l) gauge 
field. We found that the vortex flows in a direction or- 
thogonal to both the electric fleld E and the magnetic 
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FIG. 2: (color online) Phase diagram of the holographic su- 
perconductor. Bc2{T) is marked with "+"(red). Since the 
local approximation is valid below the dashed line (green), as 
shown in Ref. the Ginzburg-Landau theory is applicable 
in the shaded region. 



field i? at a constant velocity v = E/B. This is ex- 
plained by the force balance between the Lorentz force 
and the electric force observed in the conventional type 
II superconductors [T5| . 

We observed Ohmic dissipation associated with the 
vortex motion. This might be explained by the specu- 
lation that the superconducting state is violated at each 
core of the vortex lattice. In other words, the normal 
state at each core causes the energy dissipation by the 
constant motion. Since the DC-conductivity we calcu- 
lated in Sec. IV is the spatially averaged value, we can- 
not say exactly where the dissipation occurs in the vor- 
tex motion. Indeed, the dissipation is independent of the 
coefficient C(p) in Eq. (j2.13|) . It is interesting to investi- 
gate further the location of the dissipation in our model 
by calculating the DC-conductivity at each point. 

As shown in Sec. IV, the DC-current agrees with the 
current in the TDGL theory. The kinetic coefficient F 
in the TDGL equation was obtained along the line B = 
Bc2{T) in {B,T) phase diagram as shown in Fig. 2. We 
also obtained the other coefficients in the TDGL equation 
just below the critical temperature Tc at zero magnetic 
field. It is worth comparing these coefficients obtained in 
this article with the ones obtained from other phenomena 
as a consistency check. 

In general, there is a possibility that F depends on T 
and independently on B, i. c., F ~ F(T, B). To investi- 
gate the possibility, we need to evaluate F along another 
line away from the B = Bc2{T) line. It would be in- 
teresting to clarify the dependency in the phase diagram 
where the TDGL theory is available, and to compare it 
with experiments. Then, we might be able to find a sign 
of a strongly correlated condensed matter system in the 
holographic superconductor model. 
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Appendix A: vortex lattice solution 

In the zero limit of the external electric field E ^■ 
0, the solution should be reduced to the static vortex 
lattice solution obtained in Ref. Q- Let us take C(p) in 
Eq. (Pl^ as 



(Al) 



i=-c 



where pi and Ci are defined by two lattice parameters, 
ai and a2- 



Pi 



ai 



Ci := cxp [ —I — 



In the E ^ limit, 4' has a pseudoperiodicity 
*(x, y, u) = 1'(x + ai ro, y, u), 



(A2) 



(A3a) 



a2 ro 27rro 

,y , 



ai 



ai 



exp 



"27ri 


i-- 








2a J 



*(x, y, u). (A3b) 



Thus, the fundamental region Vq on the boundary 
M = is spanned by two vectors, bi — airijdx and 
b2 = a2ro/aidx — 2TTro/aidy. Here ro := 1 / V Bc2 is the 
typical length scale of the fundamental region Q . The 
triangular lattice solution, for example, is given by 



02 

ai 



oi 
2 



(A4) 



Appendix B: The vortex flow solution to TDGL 
equation 

The conventional superconductors near the critical 
temperature T — Tc arc well described by the GL the- 
ory, where Tc is the critical temperature when the applied 
magnetic field is zero. The free energy is represented by 
the order parameter O and a vector potential A as 



F= dx c{T)\DO\^ - a{T)\0\^' 



b{T) 



_D := V — ie^A, 



(Bl) 



where e, is the effective charge coupled to the vector 
potential A. Here, we assume that the parameter a{T) 
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changes from negative to positive at T = Tc as T de- 
creases, while the other parameters, b and c remain pos- 
itive near T ~ Tc- The current j is given by 



SF 



2e,c(T)(3[C'tvC']-e, lOpA) . (B2) 



The TDGL equation describing the non-equiUbrium 
states is given by a kinetic coefficient r(r) as 



SF 



= r(r) {c{T)D^ + a{T) - 6(T)|Op) O , (B3) 

where (j) is the electric potential. 

We first consider the superconducting state in the ab- 
sence of electric and magnetic fields. Setting the 1. h. s. 
of (jB3P to zero, we obtain a stationary homogeneous so- 
lution for T < Tr as 



O = Oo{T) := 



la(T) 
b{T) 



1-J ---4" m 



near the critical temperature. 

From the TDGL equation (|B3p . the perturbation 
SO := O — Oq around the homogeneous condensate (jB4[) 
satisfies the dispersion relation 



LJ = -ic{T)r{T) k 



2a{T) 



(B5) 



where we set SO oc e~*'^*+''"'^. This yields the correlation 
length ^ from the wave number fc* generating the static 
perturbation, i.e., w = as 



1 



c{T) 
2a{T) 



(B6) 



Next we consider the magnetic and electric response. 
Applying the infinitesimal magnetic field on the homoge- 
neous condensate (|B4p . the London current is generated 
by the vector potential A as [l^ 



j^-2el clT)\Oa\^A 



(B7) 



When the field strength increases beyond a critical value 
Bel, the external magnetic field begins to penetrate into 
the superconductor and vortices appear. At B — Bc2 (> 
Bel), the second order phase transition occurs and the 
superconductivity disappears. Just below the upper crit- 
ical value Bc2 , a triangular lattice appears since it is ther- 
modynamically most favorable solution (in details, see 
Ref. (3). 

For simplicity, we consider the following gauge fields 
generating the upper critical magnetic field Bc2 and a 
small electric field E in the x-direction: 



A^ -{Ba2y + Et)dx, = 0. 



(B8) 



Substituting an ansatz 

o{x, t)^ [dp c{p)e'p^ay) , y--^y + ^t (B9) 

J Dc2 

into Eq. (|B3|) . we obtain 



E 



T{T)B,2 dy 



c{T) 



r ^2 



dy2 



a(T)e (BIO) 



just below the critical temperature Tc- Here, we ne- 
glected the quadratic term with respect to ^ because T 
is very close to Tc (et is very small). Introducing a new 
variable x as 



C(y) = exp 



E 



\2c{T)T{T)Bc2 
Eq. (|B10p can be simplified as 



y x{y) , 



/ d^ _ 



aiT) 



Y := ^/2e.Bc2 



2e,c{T)Be2 
P 



(Bll) 



X + 0{E^), (B12) 



(B13) 



Neglecting the square term for small E and comparing 
Eq. (jB12|) with Eq. (j3.2bp . we finally obtain a vortex flow 
solution with the lowest energy as 



O 



dp C{p) exp ( ipx 



X exp 



e^B, 



c2 



2c{T)T{T)Be2 

2 



e*Bc 



(B14) 



with the relation 

a(T) = e^c{T)Bc2 ■ (B15) 
Substituting Eq. ((BT4)) into Eq. ((B2l) . we obtain 



], = -c{T)e,dy\OY + 
jy = c(T)e,9,|0|2 . 



2 ' ''*^ , (B16a) 



r(r) Be2 



(B16b) 



Averaging the values of the current over the a;-space, we 
find that only jx is non-zero: 



e^E 



Jx 



r(r)B, 



c2 



\o\' 



(B17) 



Eq. ((BT7| shows the Oh mic dissipation, as the electric 
field is applied in the x-direction, i.e., E ~ —dtA = Edx. 



Appendix C: The derivation of the coefficients in the 
TDGL equation 

In the following, we will determine the parameters a, 
b, and c in Eq. (|B3p just below the critical temperature 
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Tc from the correlation length and the London equation 
calculated in the holographic superconductor model Q. 

In the absence of magnetic field, the scalar field \I/ and 
the gauge potential At can be expanded as 

^{u) ^ -^{^^>^{u) + er«'2(w) + •••), 

Atiu) = a{T){q,{l -u) + eT^i{u) + 0(e|)}, (CI) 

where qc = ^L/a{Tc). Then, as shown in Ref. Q, the 
equations of motion for ^'i(u) and <E>i(u) are written by 



C4, ^ 0, 



dv? 



s{u) 



+ U + V?) ' 

(C2) 



where is defined by C in Eq. (|2.14b[) as £^ := 
—L\T=Ta- As mentioned in Sec. II, we consider the 
boundary conditions for ^1: 

lim ^i{u) = O(w^), ^i{u = 1) = regular. (C3) 

We derive the boundary conditions for from the 

requirement that the chemical potential ^ — At{Q) is 
fixed under the variation of the temperature: 

lim<li(u)=q„ $i(u = I) = 0. (C4) 

Here, the latter condition is the regularity condition 
at the horizon. The formal solution $1 satisfying the 
boundary conditions is given by 



$i(w) = — u / s{v){l — v)dv 

J u 

-{1-u) ( s{v)vdv + q^{l-u). (C5) 



In terms oi^i and $1, the correlation length ^ is rep- 
resented by 



tZ} D 



a{Tc) N 



(C6) 



in the limit T ^ Tc , where 



TV = 2 r dul ^ + $1(0) 
Jo \du ^ 

D = / du^^. 



(C7) 



Since is a solution of the linear equation (jC2[) . its 
amplitude is obtained from the next order equation: 



The boundary conditions for ^'2 are the same as the ones 
for (|C3p . Under the boundary conditions, C is Her- 
mitian for the inner product (|4.10p . This yields 



-2qc{ 



1 + U + 



(C9) 



Substituting Eq. (|C5|) into Eq. (jC9[) . we obtain the ampli- 
tude A defined by A := for a normalized solution 
■01 satisfying V'i(l) = 1 as 



A' 



S2(0) 



S„(u) 



lo mu)du - i^J^ Si(w)du 
1 (I - 



-dv 



= 2). (CIO) 



Numerical calculation determines the value of the am- 
plitude A and hence C2 in Eq. (|2.5p is evaluated as 



C2 



lim ^ 



6.55 e. 



1/2 



^ in Eq. (|C6p is also evaluated as 



-1/2 



i ~ 0.0635 X 



T, 



(Cll) 



(C12) 



For the vector potential A generating small magnetic 
field, the London equation just below T = Tc Q is eval- 
uated as 



lim ( J 



du 



IIP 



2 hi^ \ Jq C9 



c^A 



X 1.172 X 10" 



O2 



(C13) 



where we used Eqs. (|2.1bp and (|2.6a[) to derive the third 
equality. Comparing Eq. (jC13[) with Eq. (|B7p and iden- 
tifying ( O2 ) with Oq, we obtain the coefficient c in the 
TDGL equation (|B3P just below Tc as 



2nl 5.85x10"'' 



L2 



T3 



(C14) 



where we used the fact, e, = 1 derived in Sec. IV. Substi- 
tution of Eqs. ((CT2I) and ((CT4)) into Eq. (|B6]) yields the 
coefficient a just below Tc as 



^^2 



1 + u + 



(C8) 



a(T ^ T,: 



2k| 0.0726 ct 



(C15) 



^ Under the asymptotic boundary condition in Eq. ||C4| |. e 
q/qc — 1 defined in Ref. Q is equal to ey. 



h is also evaluated from Eqs. (|2.6a|) . (|B4|) . and (|C11|) as 
-2 \ 3 o 7c; V 1 n-6 



K^c) 



/2«4 



X 



2.75 X 10- 



(C16) 
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Appendix D: Ohmic dissipation 



term in the above equation as the Ohmic dissipation: 



Since the bulk spacetime possesses Kihing vector = 
= (9t,)^, the energy-momentum tensor of the form 



- ^ 



9^u 



^2 



(Dl) 



satisfies the conservation law, V^(T'^j,^'') — 0. Thus, 
we obtain 



0=. j d^x^gV^ (T^r) 



+ / ds^T^r 

H 



(D2) 



bdy 



where, '^f,i, 'H, and bdy represent v = const, null hyper- 
surfaces, the null hypersurface at the black hole horizon, 
and the timelike AdS boundary, respectively. Since the 
spatially averaged value of T^i, does not depend on u, the 
first and the second terms in the second line of Eq. (jD2[) 
cancel each other. This implies 



d^aT\e 



H 



d^aT\e 



(D3) 



bdy 



Due to the rapid fall-off condition for p.l6p . only 
the U{1) bulk gauge field contributes to the boundary 



/ dY.aT°-U^ = / dvd^x{ J' ) E, 

Jhdy Jhdy 



'bdy 

Here, we used the fact 



-9 



u=0 



(i = x,y) 



Hence, Eq. (jD3|) is reduced to 



2„2 , 
4^ JH 



dv d^x S'^ F„Fyj +2{L ay 



bdy 



(D4) 



We can extract the subleading terms at 0(e) from 
above. Noting 



dxdy djay = , 



dxdy fyj — / dxdy 2d[yaj] 



we obtain 



dvd'^x Dy^ 



n 



bdy 



dv d^x S{J') E, 



(D5) 
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